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Abstract 

We discuss constraints on three flavor neutrino mixings from the accel- 
erator and reactor experiments, the Kamiokande multi-GeV data, and the 
solar neutrino observations. The LSND result is excluded at 90%CL by the 
constraints imposed by all the data of reactor and accelerator experiments 
and the Kamiokande multi-GeV data if the mass scale required for the solu- 
tion to the solar neutrino problem is hierarchically small. The region of a set 
of the effective two-flavor mixing parameters (Am 2 , sin 2 29) is given for the 
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channel — > u e which is allowed at 90%CL by the multi-GeV Kamiokande 
data alone. 
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Recently the LSND collaboration has claimed that they have found candi- 
date events for — > v e oscillation (See also [0]). If their result turns out to 
be correct, it gives us an important information for masses and mixing angles 
of neutrinos. In this paper a possibility is explored that all the experimental 
data, solar, atmospheric, reactor and accelerator data including LSND can 
be explained within a framework of three flavor neutrino mixing. It turns out 
that the LSND result is excluded at 90%CL by the constraint imposed by all 
the data of reactor and accelerator experiments and the Kamiokande multi- 
GeV data. The statement is true even without invoking actual solar neutrino 
data as far as the mass scale required for the solution is hierarchically small. 

Among the atmospheric neutrino data § § @ [| we discuss 
only the Kamiokande multi-GeV data |4| throughout this paper. The rea- 
son for this restriction is three-fold; (1) It is the only experiment that gives 
us a nontrivial zenith-angle dependence which leads to the upper and lower 
bounds for the mass-squared difference of neutrinos. (2) The result of Monte- 
Carlo simulation for neutrino energy spectrum is published only for the the 
Kamiokande multi-GeV data. (3) It seems difficult to reconcile it with NU- 
SEX § and Frejus % data. 

The Dirac equation for three flavors of neutrinos in vacuum is given by 

= UdmgiE,, E 2 , E 3 )U-^(x), (1) 

which is easily solved: 

= E^iC^**^), (2) 

where Ej = ^Jp 2 + mj (j = 1,2,3) is the energy of neutrinos in the mass 
basis, ty a (x) = (v e (x), ^(x), z/ T (x)) T (a—e,/i, r) is the wave function of 
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neutrinos in the flavor basis, and 



U = 



( U el U e2 U e3 
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with Qj = cos 9ij, Sij = sin^j is the orthogonal mixing matrix of neutrinos. 
We will not discuss the CP violating phase of the mixing matrix here for 
simplicity. The probability of u a — > vp transition is given by 



P{v a -> vp\ E, L) 



4 sin 



2 / \ V 2 V 2 
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for a = (3 
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where AE^ = E t - E j ~ (mf - m 2 j )/2E = Am 2 ij /2E is the difference of the 
energy of two mass eigenstates. 

The number of neutrinos vp {{3 = e, /i, r) is measured in terms of charged 
leptons tp which come out from a scattering vpX — > ipX' . In appearance 
experiments of v a — > vp the expected number of the charged leptons ip is 
given by 

N{y a ^tp-L) 

= n T r dE f maX dq e(q)E a (E) dap( f' q) P(v a -> up; E, L), (5) 
Jo Jo dq 

whereas in disappearance experiments we measure attenuation of beam 

N aa (L)-N(v a ->£ a ,L) 

= n T r dE T maX dq e(q)F a (E) dCTa{ f' q) (1 - P(u a -> i/ a ; E, L)) , (6) 
jo Jo dq 
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where 

N a(3 (L)^n T J Q dEj Q dqe(q)F a (E) ' q) . (7) 

F a (E) is the flux of neutrino v a with energy E, tit is the number of tar- 
get nucleons, e(q) is the detection efficiency function for charged leptons £p 
of energy q, da p{E , q) / dq is the differential cross section of the interaction 
vpX -> ipX', P(v a -> i^; £?) = |(z/ a (0)|z/ /3 (L))| 2 is the probability of u a -> ^ 
transitions with energy E after traveling a distance L. The results of exper- 
iments are usually expressed in terms of a set of the oscillation parameters 
(Am 2 , sin 2 29) in the two-flavor analysis. The probability in the two-flavor 
mixing is given by 



P{y a -> v fi \ E, L) 

1 - sin 2 29 a sin 2 ^ ~^r~ ^ for a = /5 

2 „ 2 / Am 2 L\ 
sm 26» Q!/3 sin [ 4 ^ ) for a ^ (3. 



(8) 



Introducing the notation 
'Am 2 L\\ 



sin 2 



4E 



_i „ p / \ f (F , d<T P (E,q) 2 ( Am 2 V 



N ( 

we have the expected number of charged leptons in the two-flavor scenario 
, NW(v„^e a ;L) . 2 „. ,. 2 J . 2 ( Am 2 L \\ 

1 JUST- = sm * (Am 1 \ sm l^HL. 

(disappearance experiments) 



(9) 



NQ\v a -^i p ;L) . 2 2 / , 2 / Am 2 V 
— = sm 29 a ,(Am ) (s.n 



(appearance experiments), (10) 
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where we have denoted the Am 2 dependence of the mixing angle 8 (Am 2 ) 
explicitly to indicate that 8(Am 2 ) is a function of Am 2 in the two-flavor 
analysis. 

The boundary of the excluded region (for negative results) in the (Am 2 , sin 2 26) 
plot is determined by 

NQ\v a ^l s ;L) . 2 2 J , 2 ( Arn 3 L \ \ 

f = N all (L) = S1 " 2 ^ (Am » \ S1 ° [IE-) L e ' (U) 

for appearance experiments where the charged leptons are detected at one 
point at a distance L, or 

# (2) (f«-*«;£i) N^(u a ^£ a ,L 2 ) 



NaaiU) N aa {L 



2 j 



sin 2 26 Q (Am 2 ^ 



;i2) 



for disappearance experiments where the charged leptons are detected at 
two points at distances Li and L 2 (Li < L 2 ). In (JTTj) and (|T^) e denotes 
the largest fraction of the appearance events allowed by a given confidence 
level, i.e., N^^a — > £p;L)/N a p(L) < e for appearance experiments, and 
the largest fraction of beam attenuation, \N^ 2 \u a — > £ a ; Li)/N aa (Li) — 
(i/ Q — > L 2 )/A^ QQ ,(L 2 )| < e for disappearance experiments. 
From (|TT|) and (0), we can read off the value of (sin 2 (Am 2 L / AE)"j for ar- 
bitrary Am 2 from the figure of the two-flavor mixing parameters (Am 2 , sin 2 29) 
given in each experimental paper as long as sin 2 26(Am 2 ) < l 3 : 

Am 2 L 2 \\ /. 2 (Am 2 Li\\ e 



sin 2 — — — I > — ( sin 2 



AE \ \ AE sm 2 26JAm 2 ) 

(disappearance experiments) 



3 In case of the CP violating phase, which will not be discussed in the present paper, 
this is not the case. To discuss the CP violating effect, one needs the information of 
^sin ( A 2 e £ ) ) ' wmc h cannot be obtained from the information in published papers only. 
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sin 2 



/ Am 2 // 



V AE J/ a ^ sin 2 2# Q/3 (Am 2 ) 

(appearance experiments). (13) 

In ([13]) and in the following, sin 2 2# a (Am 2 ) or sin 2 2^(Am 2 ) stands for 



the value of sin 2 2#(Am 2 ) on the boundary of the allowed region in the 
(Am 2 , sin 2 Iff) plot in the two-flavor analysis. 

Similarly we can express the number of the expected charged leptons in 
the three flavor mixing: 



From C|l3f ) we observe that the quantity (sin 2 (Am 2 L / AE)^ /e is equal to 
sin 2 2$ (Am 2 ) which can be read off from the published literatures, and we 
can express the conditions for the three flavor mixing parameters in case of 
negative results: 

iV (3) K^4;£i) _ iV (3) K->l a ;L 2 ) = v e 2 2 

N aa (Li) N aa (L 2 ) ^sin 2 20 Q (Am?.) ai ^ 

(disappearance experiments) 
N&(v a ^l p] L) „ e rr rr rr rr 

£ > tf^^'^^ 

(appearance experiments). (15) 

Notice that the left-hand side of fll3|) is defined independent of the number 
of flavors of neutrinos. 

Throughout this paper we assume that a single mass scale is involved 
in the solution of the solar neutrino problem, which is hierarchically small 
compared to others. Namely, we assume that 

Am 21 < Am\ 2 < Am^, (16) 
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where we have assumed m\ < m\ < mj without loss of generality. In the 
case where Am 32 <C Am^ < Am^, we can show that we obtain the same 
conclusions, although we will not give the calculation here. 



The hierarchy ( |l6l) is satisfied in the two-flavor mixing solution to the 
solar neutrino problem || |10| which requires, 



(Am 2 , sin 2 29) 
(Am 2 x , sin 2 29 12 ) Q 

(O(10- 5 eV 2 ),O(10- 2 )), 
(O(10- 5 eV 2 ),O(l)) 



(vacuum solution) 
(small angle MSW solution) 
(large angle MSW solution). 



(17) 



These mass scales are much smaller than those which appear in the atmo- 
spheric neutrino observations ||, or in the LSND experiment 0. In fact the 
only condition on (Am| 1; sin 2 29\ 2 ) which is essential to the discussions below 
is (|T6D, and our conclusions remain unchanged irrespective of the value of 812 
as long as ( pf ) holds. 

One can show |12|] [13|] that under the mass hierarchy ( |T6| ) the relation 



U 2 e3 « 1, 



(18) 



must hold in order to have solar neutrino deficit under the constraints from 
the accelerator and the reactor experiments. In this setting it can also be 
demonstrated that the solar neutrino problem is indeed solved by a two-flavor 



framework in the MSW and the vacuum solutions 12 



In the present case the formulas in (14) become much simpler [11]: 

u a ) = 4 sin 2 ( ^] Ul x Ul 2 + 4 sin 2 ( ^] U^l - t/ 2 3 ) 



1-P(v n 



P->, - up) = -4 sin 2 ( U al U a2 U, x U, 2 + 4 sin 2 ( UlM 



AE 



AE 



(19) 
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For the cases discussed below, we have \Am\ 1 L / AE\ <C 1 and the first term 
in each formula in ( PUD can be ignored. Hence we get 

4[/^ 3 (l — U^g) < sin 2 2# a (Am 31 ) (disappearance experiments) 

4C/ 2 3 C/| 3 < sin 2 26 , Cfj g(Am 31 ) (appearance experiments) (20) 

for negative results, respectively. 

The data by the LSND group suggests that the allowed region for the 
mass-squared difference Am| x is approximately Am| x > 5 x 10~ 2 eV 2 . If we 
combine the data on the same channel — > V e (and — > u e ) by the E776 
group [14], however, the region Am 31 >2.5eV 2 seems to be almost excluded, 



and we have 5 x 10 2 eV 2 < Am^, < 2.5eV 2 at 90% confidence level. 

As has been pointed out in Ref . [|ll|] , strong constraints on the mixing an- 
gle come from the reactor experiment ||15|| . Using (|20"D we have the constraint 
from the reactor experiment Jl5| 



sin 2 20 13 = 4t/ 2 3 (l - U 2 e3 ) < sin 2 2# B uge y (Am 2 1 ), (21) 

where sin 2 2#Buge y (Am 2 a ) stands for the value of sin 2 26 on the boundary of 
the allowed region in the (Am 2 , sin 2 26) plot in lfL5]| . For the entire region 
5 x 10~ 2 eV 2 < Am 2 1 <2.5eV 2 , sin 2 2fl B ugey(Am 2 1 ) is small pf: 

^<sin 2 2^Bu g c y (Am 2 1 )<^. (22) 

From (^l|) we have either small t/ 2 3 or large £7 2 3 , but the latter possibility is 
excluded by ([Tj^). Therefore we have 

4s = U 2 eZ < \ sin 2 2# Bugey (Am^). (23) 

On the other hand, we have another constraint from the disappearance 
experiment of |16] 



4s 2 3 c 2 3 (l - S 2 3 c 2 3 ) = 40* (1 - U%) < sin 2 2^ C DHSw(Am 2 1 ), (24) 
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where sin 2 26 , cDHSw(Am 2 a ) stands for the value of sin 2 29 on the boundary 
of the allowed region in the (Am 2 , sin 2 29) plot in |]IE||. The mixing angle in 



this case is constrained for 0.7eV 2 < Ami < 13eV 2 as 

sin 2 2^ CDH sw(Am 2 1 )<0.2. (25) 

If we consider the probability Piv^ — > v^) in the atmospheric neutrino ex- 
periments for the mass region above, we have small deviation of P[y^ — > v^) 
from unity 

1 - PK - ^) * 2t/ 2 3 (l - U^)<0.1, (26) 

where we have averaged over rapid oscillations. (|26|) is obviously incon- 
sistent with the atmospheric neutrino observations || ||] § [§], since we 
cannot have gross deficit of in this case. In fact we have verified explic- 



itly, by the same calculation as in Ref. |T7|], that the region Am^ > 0.47eV 2 



with the constraints (pl"|) and ( f2~4"|) is excluded at 95% confidence level by 
the Kamiokande multi-GeV data (x 2 >17.6 for 3 degrees of freedom, where 
X 2 lin =3-2, (x 2 ~ xLn)/7 - 2-1 implies 2a). 

So we are left with the region 5 x lCT 2 eV 2 < Am| x < 0.47eV 2 . Applying 
the formula ( |T9"D to the case of LSND, we have 

sj 3 sin 2 29 13 = 4U 2 e3 U 2 ^ = sin 2 20 LSND (Am^), (27) 

where sin 2 26 , LSN£)(Am 2 1 ) stands for the value of sin 2 29 within the allowed 
region in the (Am 2 , sin 2 29) plot in and the LSND data indicates 

1.5 x 1(T 3 < sin 2 2^ LSND (Am 2 1 ) < 1. (28) 

From (^) and the constraint s| 3 < 1, it follows that 

sin 2 2# LS ND(Am 2 1 ) < sin 2 2^ B u g e y (Am 2 1 ), (29) 



10 



so that we have 



Am 2 31 >0.25eV 2 , 



(30) 



with 



1.5 x 1(T 3 < sin 2 2^ LSND (A 2 1 )<3.8 x 10" 2 , (31) 



where we have used the (Am 2 , sin 2 26) plots in [15] and [l|]. However, we have 



verified explicitly again that the region Am 2 > 0.25eV 2 is excluded at 90% 
confidence level by the Kamiokande multi-GeV data (x 2 >15 for three degrees 
of freedom, (x 2 — Xmm)/7 — implies 1.6a). Therefore, we conclude that 
the LSND data cannot be explained by neutrino oscillations among three 
flavors, if all the accelerator and reactor data as well as the Kamiokande 
multi-GeV dat db £1X6 taken for granted. 

The present LSND data allows conflicting interpretations either as a pos- 
sible evidence for neutrino oscillation [|I|] , or a stringent bound for the mixing 
parameters . It might be possible that the allowed region of the set of the 
parameters (Am 2 , sin 2 26 l ) implied by the LSND data changes in the future. 
We have looked for the region of (Am 2 , sin 2 26) for — > v e (or 77 M — > ]7 e ), 
in which —>■ v e oscillation is consistent with all the experiments (except 
LSND), including the Kamiokande multi-GeV data. To obtain the (Am 2 , 
sin 2 26 plot of the two-flavor analysis for general — > v e oscillations with 
a mixing angle ^(Am^), we use the following correspondence between the 
rates for the v„ — > e process in the three and the two-flavor frameworks: 




(32) 



11 



where we have used the hierarchical condition fllED- From and fl32|) we 
have 



sin 2 26^ e (Am^) ~ s 2 3 sin 2 2#i 3 

< S23sin 2 26»Bug e y(Am^). (33) 

The Kamiokande multi-GeV data has better fit for larger values of sin 2 2# 23 , 
so (|3~3"|) shows that sin 2 26 l Aie (Am| 1 ) cannot be larger than sin 2 26 , B U ge y (Am 2 1 )/2, 



if sin 2 2# 13 <C if]. The result is shown in Fig.l. 
(Insert Fig.l here.) 

The region suggested by the LSND experiment JE] is close to the 90%CL 
region obtained in our analysis. If the 20% systematic uncertainty mentioned 
in Ref. U shifts the allowed region of LSND in the direction of smaller 
mixing, there may be a chance that all the neutrino anomalies are explained 
by three flavor neutrino oscillations. Hopefully further data from the LSND 
group will clarify the situation. 
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4 This is not the case for Am 2 ^ <3 x 10 2 cV 2 , where 1 — — > v^) = 4s 2 , 3 c 2 3 (l 



S23CJ3) ~ 0(1) becomes possible even for 623 < 7r/4. 
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Figures 



Fig.l Below the solid and dashed lines are the regions for — > v e (or 
— > 17 e ) oscillation which are compatible with all the experiments (ex- 
cept LSND), including the atmospheric multi-GeV data of Kamiokande 
at 68%CL (solid) and 90%CL (dashed), respectively. To get this plot 
of the two-flavor mixing parameters, we have used the correspondence 
sin 2 29 = sin 2 2^ e (Am 2 1 ) = s 2 3 sin 2 2df 3 (cf. flU). The shadowed area 
stands for the region allowed by all the accelerator and reactor experi- 
ments including LSND (The dotted, dashed and dot-dashed lines stand 
for the LSND E776 [14]], and Bugey |1| experiments, respectively). 
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Fig.1 



